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Abstract 

We suggest a duality between the statistical and standard (dynamical) dis- 
tributions of partons in the nucleons. The temperature parameter entering 
into the statistical form for the quark distributions is estimated. It is found 
that this effective temperature is practically the same for the dependence on 
longitudinal and transverse momenta and, in turn, it is close to the freeze-out 
temperature in high energy heavy-ion collisions. 
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1. Introduction 

As is known, the statistical (thermal) models have successfully been ap- 
plied to describe ratios of hadronic yields produced in heavy-ion collisions 
(see, for example, [2I, and references therein). At the same time, 

the source of very fast thermalization is currently unknown and alternative 
or complementary possibilities to explain the thermal spectra are of much 
interest. 

Situations where statistical models were applied while the notion of sta- 
tistical system was not obvious are not uncommon. In particular, the statis- 
tical model was successfully applied to analyze deep inelastic lepton-nucleon 
scattering (DIS) 0, 0, where the statistical form j^, [oj for distributions of 
unpolarized and polarized quarks and gluons (partons) in a nucleon was used. 

This model may be compared with more standard parametrization of 
parton distributions based on the Regge theory at low x and quark counting 



rules at large x [10|]-[l2|. It is especially interesting to find a counterpart 
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and understand a physical meaning of the universal temperature introduced 
in the statistical model. This is the main goal of this paper. Making this 
comparison we suggest a concept of the duality between the statistical and 
dynamical descriptions of the DIS and strong interactions of particles and 
explore its consequences and implications. 



2. Statistical and dynamical description of quark distributions 

2.1. Quark distribution in proton: longitudinal momentum 

According to many experimental data, inclusive spectra of hadrons pro- 
duced in heavy-ion collisions can be fitted using a thermodynamic statistical 
distribution as a starting point for the system of final hadrons, see for exam- 

pie [il,8 

= C^{exp((e.-/..)/T) ± ir\ (1) 

where -|- is for fermions and — is for bosons, et and fit are the kinetic energy 
and the chemical potential of the hadron h, T is the temperature, and is 
the normalization coefficient. For mesons simplifying this case we can assume 
that /i/i ^ 0, (in fact, it generally cannot be strictly zero [^). At high energies, 
when e/i/T >> 1, Eq.([T]) is usually presented in the Boltzmann statistic form 

^ C^exp(-e,/T) . (2) 

A statistical approach was applied in great detail to study the quark struc- 
ture of nucleons in Ref. 0, using earlier approaches as a starting point js], [oj. 
The nucleon is viewed as a gas of massless partons (quarks, antiquarks, glu- 
ons) in equilibrium at a given temperature in a finite size volume. According 
to the thermodynamical approach, the quark distribution in a proton q{x) 
in the infinite momentum frame (IMF) at the input energy scale Qq, which 
is used in the perturbative QCD [l3| as the initial condition for the evo- 
lution, is fitted as a function of the light cone variable x in the following 
form: 

""^^"^^ exp {{x - Xog)/x) + 1 exp(x/x) + 1 ' ^ ^ 

where x plays the role of a temperature inside the proton and Xoq is the chem- 
ical potential of the quark inside the proton. All the parameters are found 
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from the description of the DIS. The same form is suggested in 0] for anti- 
quarks xq{x) by changing the sign for the chemical potential Xq^. Therefore, 
the proton structure function ^2(0;, Q^) = Yl, flavour-^ i^lfi^^'Q'^) + Qf{^iQ'^)) 
has also a thermodynamical form like Eq.Q. 

Note that the statistical weight of the quarks can be written in the form 
exp{{Eq — /i)/T), where Eg = {P ■ Pq)/m is a quark energy in the nucleon 
rest frame, Pg and P are the four-momenta of quark and nucleon, respec- 
tively, m is the nucleon mass. Passing to the infinite momentum frame the 
quark distribution in a nucleon over the longitudinal momentum pz can be 
represented in the following statistical form [9|: 

where E, P and m are the energy, the momentum and the mass of the 
nucleon moving in the z-direction, e and Pz are the energy and momentum 
in the ^-direction of a parton in the nucleon, T is the temperature and yU is 
the chemical potential. For massless partons and /i ~ one gets 

q{x) ~ exp (-^) = ^^P(~|) • (^) 

where x = Pz/P is the longitudinal fraction of the parton momentum and 
X = 2T/m. This form agrees with the result obtained in js], [oj. The inclusion 
of of the quark transverse momentum kqt results in Ee — Ppz — m'^x{l + 
^qt/ {x'^fTT''^))- It means that the relativistic invariant variable x is replaced by 
the variable x{l + kg^/{x'^m'^)) Q. 

The statistical parametrization of parton distributions [gI, 0] may be com- 
pared with a more standard form (see, e.g., [loj) describing the experimental 
data on the DIS: 

xq{x) = AsX-'^'il - xf' + B^x^'^il - xf" , (6) 

where the first term is due to the sea quark distribution which is enhanced 
at low X, according to the DIS experimental data, whereas the second term 
is due to the valence quark distribution. Note that to fit the experimental 



data on the DIS a more complicated form for xq{x) is also used [ll|, |12 . 

In particular, the valence part ai x < l/bg^ can, to a very good approxi- 
mation, be written as: 

xq^{x) ~ Bg^x'^" exp {-xbg^) . (7) 
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Figure 1: The valence u-quark distribution in the proton as a function of the Bjorken 
variable x at Qq — 4:{GeV/cy . The solid curve is the CKMT parametrization [lO|, the 
dashed line corresponds to the statistical parametrization given by Eq. Q the dotted curve 
corresponds to the BS statistical (thermal) model suggested in pil3|- 



This form is similar to Eq.(l5]) ii x = In this case the proton tempera- 

ture is T = X ■ m/2 = m/{2hv{Q'^)) for current quarks which are practically 
massless. Approximately one has for the valence w-quarks h, ~ 2.5 — 3 
and for d-quarks ~ 3.5 - 4 at = 2 - 4 (GeV/c)^ Gj; therefore, 
T ~ 120 — 150 MeV for the massless quarks. 

The similarity of the thermal form for the quark distribution given by Eq. (|5]) 
and the dynamic form for (Eq.(I7]) maybe interpreted as a duality of the 
thermal and dynamical descriptions of the parton distribution in proton. This 
property is by no means surprising as both parametrization are fitted to de- 
scribe the same data. In Fig.l, the valence w-quark distribution in the proton 
as a function of x is presented, the solid curve corresponds to the CKMT 
parametrization [lo[, see the second term of Eq.([6]) at Qq = AiGeV / cY] the 
dotted line corresponds to the statistical (thermal) BS model [TJ , see the first 
term of Eq.([3]) at the same value of Ql; the dash-dotted curve corresponds to 
the statistic parametrization given by Eq.(I7]). Figure 1 shows that all three 
lines are very close to each other at 0.01 < x < 0.4. 

Let us compare the parameters of the statistical parametrization for u.u{x) 
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given by Eq.([7]) and the BS parametrization (Eq.([3])). While the temperature 
has an evident dynamical counterpart T = m/{2by) , a similar relation for 
the chemical potential Xoq of valence quarks is not obvious. The temperature 
T = m/{2hy) entering into Eq.Q is about 120 — 150 MeV, whereas the model 
of Ref. [71] suggests that T is about 50 MeV and the chemical potential for 
u-quarks Xqu — 216 MeV as first found in [sl, 0]. 

2.2. Quark distribution in proton: transverse momentum 

The quark distribution function in a proton using the variable x and the 
transverse momentum kt can be represented in the factorized form fq{x, kt) = 



q{x)g{kt) which of course cannot be valid at all values of x 15|. This is sup- 



ported by simulations within the lattice QCD [16| and by the observation of 



the so-called "seagull" effect [17[, i.e., the weak x-dependence of the mean 
transverse momentum < > of hadrons produced in hadron inclusive reac- 
tions at low X, e.g., X < 0.5 and the strong x-dependence of < > (x) at 



X — )■ 1, see for example [18j and references therein. 

Developing the approach of the previous section one can fit the quark 
distribution g{kt) also in the statistical form ([2]), i.e., 

g{h) ~ exp(-|i) , (8) 



where e^j = yk"^ + is the transverse energy of quarks in proton and rUg 

is the quark mass. For the massless quarks (the current quarks) g{kt) can be 
represented in the following form: 

kt kt 
g{kt) ^ exp{-—) = exp{-—^). (9) 

Applying here the same effective temperature T ~ 120 — 150 MeV, as for 
longitudinal momentum distribution, we get similar results on < kt >, which 



were obtained for valence current quarks (see e.g. [IJ] and Ref. therein). 
Actually, these values for < kt > of quarks in proton are used to get the 
inclusive transverse momentum spectra of hadrons produced in p—p collisions 



at not large pt, which were obtained in 18|]-[22[ within the dual parton model 



(DPM) [Ig or the quark-gluon string model (QGSM) |20l. 

Therefore, the transverse momentum distribution of partons in proton 
can also be described in the statistical form with the same value of the 
temperature T as the parton distributions in proton over the longitudinal 
momentum or over its fraction x. 
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Note that the longitudinal and transverse momentum dependences have 
been related in a model approach implying Lorentz (rotational) invariance 
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3. Inclusive hadron spectra in N — N and A — A collisions 

The inclusive spectra of hadrons produced in — collisions at high 
energies in the central region are not described by the statistical Fermi-Dirac 



distribution given by Eq.(IT]) (see, e.g., the fit of experimental data in [24 
and references therein). For example, the transverse momentum pt spectra 
of all the hadrons produced in p — p collisions at high energies have a more 
complicated form than Eq.(IT]). However, the experimental data on these 
spectra at not large values of pt can be described [l8]-[22|| within the QGSM 
using an exponential form similar to Eq.([9]) for the quark distribution in 
proton on the internal transverse momentum fcj. 

This manifests another application of the duality principle suggested in 
Section lb. 

Contrary to this, the experimental ratios of hadron yields and trans- 
verse mass spectra of hadrons produced in central heavy-ion collisions are 
described rather satisfactorily by the statistical (thermal) model, being its 
traditional field of application (see, for example [H, 25, 0], j26[- 28|, and ref 



erences therein). Note that the freeze-out temperature parameter obtained 
in this analysis at the zero chemical potential is also compatible with our 
estimation T ~ 150 MeV. 



4. Conclusion 

We analyzed the parton distribution functions in a proton using the sta- 
tistical model jsl, [s], lif. We suggest a duality principle which means a sim- 
ilarity of the thermal distributions of partons (quarks and gluons) and the 
dynamical description of these partons. This duality allowed us to find an 
effective temperature T ~ 120 — 150 MeV for the massless quarks. These 
values for T are the same for the distributions of partons on the longitudinal 
and transverse momentum. 

The inclusive spectra of hadrons produced in A^A^ collisions at high ener- 
gies are not described by the thermal statistical distribution. However, using 
the thermal statistical form for quark distribution as a function of the longi- 
tudinal momentum fraction x and the internal transverse momentum kf one 
can get a satisfactory description of the experimental data on these spectra. 
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It happens that the freeze-out temperature in central heavy-ion colhsions 
at zero chemical potential has a similar value which was estimated in this 
paper for the massless quarks in a free nucleon. 
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